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THEORETICAL STUDIES OF SOI^-PUMPED LASERS 


By 

Wynford L. Harries* 

SUMMARY 

Estimates of the absorption and emission characteristics of molecules 
are required to predict new materials for solar-pumped lasers. These char- 
acteristics can be described terms of the Franck Condon factors, idiich 
are calculated from the molecular wave functions. 

Wave functions for vibrational levels in the lower and upper electronic 
states of I 2 and NSe are calculated numerically and methods of checkinpc 
errors discussed. Errors arise tdien the vibrational quantum numbers are 
high; but, using a calculated rather than measured value of the dissociation 
energy, wave functions up to the fiftieth vibrational level are obtained. A 
new numerical method of evaluating the wave functions is given, which should 
be more accurate in the region of electronic transitions during absorption. 
Franck Condon factors, plotted versus the wavelength of the absorbed 
photons, are shotm, and a check on the Franck Condon factors is made using 
the vibrational sum rule. 


INTRODUCTION 

In order to predict the performance of solar-pumped lasers as power 
converters, a knowledge of the absorption and emission characteristics of 
molecular materials is required. The data on several thousand diatomic 
molecules have been tabulated by Herzberg (ref. 1) and Huber and Herzberg 
(ref. 2). The information, where complete, yields the potential energy 
curves and vibrational levels of the molecules. The molecular wave func- 
tions can then be calculated and the Franck Condon factors obtained, which 
give information on the absorption and emission characteristics. The 

♦Professor, Department of Physics, Old Dominion University, Norfolk, 
Virginia 23508-0369 


c«lcul«tion of the wave functions is outlined in the next section, "Molec- 
ular Functions," which includes a discussion of the errors and tsethods of 
checking them. A new numeric method of calculating the wave functions 
using an asymptotic solution is also presented. There are many advantages 
to this new approach, including greater accuracy in calculating the 
probability of electronic transitions. Under the subsection titled "Franck 
Condon Factors," the Franck Condon factors are obtained from the wave func> 
tions and are calculated for a given lower and upper level. The energy 
difference between lovrer and upper levels corresponds to the wavelength ^ 
of Che absorbed photon; for a given lower and a series of upper levels, the 
Franck Condon factors are calculated and plotted vs. X. The resulting 
curve is an approximate indication of Che absorption cross section (per unit 
wavelength) vs. X. The final sections of the text discuss conclusions and 
work to be performed. 


MOLECULAR WAVE FUNCTIONS 
Introduction 

The molecular wave functions were first calculated following Che method 
of Fazio (ref. 3). The wave functions are those of one atom of a biatomic 
molecular in a frame of reference whose origin is centered on the other 
atom. They are obtained as solutions of Schrodinger's equation in the 
radial form. The atom in question is assumed to be in a potential well 
V(r) given by a Morse potential energy function (ref. 4): 


V(r) ■ D (1 - exp t-6(r - ri)}^ 
e 


( 1 ) 


where Dg is the dissociation energy, r^ the equilibrium internuclear 
distance (fig. 1), and 0 is given by 


2 


0 - ( 


2 c M 
h D 


.) « 


Here c is the velocity of liRht, y is the reduced ness of the stom, h 
is the Planck constant, and is the harmonic molecular frequency. 

Substituting equation (1) into the Schrodinger equation yields 


.UlijL [e - 
k2 L 


dr^ h 


D [l - exp (-0 ( 

ft 


r - r,))]^] 


♦ - 0 (3) 


The energy E is the eigenvalue for the anharmonic oscillator (ref. 


E(v) « hcw^ [(v ♦ i) - (v ♦ (v 


where is in cm"^ (cUg “ v, the frequency) and v is the vibra- 

tional quantum number. The first term is the solution for a harmonic oscil* 
lator, the last two arise because of the anharmonic nature of the Mbrse 
function. 

The quantities u, r^ , WgXg, (»e^e available (refs. I, 

2) for many materials. For the lower level D^, the energy difference 
between the ground vibrational state and the dissociation energy is given 
(fig. 1); hence, 


D ■ D ♦ -r hew 
e o 2 e 


Th.? value of w^y^ is sometimes included in reference 1 but not in 
reference 2. 
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For Che upper electroaic level, there exists enocher set of constants 

I 

^e*' '*’e^e'* ^eFe* * 'l* mininua energy of 

Che upper curve T^, but D2 for Che upper level (fig. 1) Is calculated 

from 


( w ^?)2 

■ 4 al' X * 
e e 


(5) 


which assumes Chat o>«yze' * 0. Usually cj^ » «>eXe *^6^6 ^e' 

“eXe' » ^eFe* values of equation (4) is accurate. For 

For V > 10 appreciable errors occur because Che cubic term then becomes 
important . 

The wave functions are obtained ntimerically (refs. 4-6). The value of 
E is given by equation (4) for a chosen v. The classical turning points 
of the potential well are then obtained from equation (1): 




ri 




( 6 ) 


For r » R^, Che coefficient of 4' in equation (3) approaches a constant, 
and there is an asymptotic solution: 


♦ ♦ exp 




-2 


5.951 X 10 


U [v(r) - E)] 



(7) 


and here 


dr (r - ri ) 


( 8 ) 
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SubtciCution into the Schrodinger equation yields: 


^ - 5.951 * lo“ u (v(r) - b) (tj. ♦ 4^ |I) 

dr2 " 2 


where dr is an increment towards the origin and negative. Values of 
and (d^/dr) at a distance dr closer to the origin are 


dt|> d^'P dr 


(d*j 

'•dr ^2 dr 


The scheme proceeds iteratively to smaller r, and, as noted previ' 
ously (ref. 4), the function diverges when r < R_. By monitoring the 
signs of il* and d^Jdr, the computer recorded * 0 tdien divergence 


occurred . 


The normalized function obeys 


!Z ♦^(r) dr - 1 . 


The values of t|i were stored, and ^ integrated by Simpson's rule. Values 


were then obtained of the normalized function 


Vr) 


f/® 'l^(r) dr 


5 




whtre Che llmlcs A end B ere choeen to saclsfy en accuracy of five 
significant figures* The plots of tN superloiposed onto tho Morse 
curves (fig* 2)* The poelclon of the wave functions on the ordlnece 
corresponda to the energy of the level, but their amplitudes are arbitrary 
although relative to each other* The reiteration scheme started for the 
lower curve at an upper limit of r such ths.t t < 10"^^ and proceeded to 
a lower limit which was less chan R_ in 2,490 increments* Every tench 
value was stored for integration* For the lower level, the scheme started 
at r ■ 3*318 A and proceeded to 2*318 A* For Che upper level, the 
integration was carried out in two parts: for the region from 5*051 to 

3*318 A in 2,490 steps and from 3*318 to 2*318 a in 2,490 steps* The two 
regions were required for normalisation, out Che limits of r in the latter 
regions coincided with those of Che lower level so that the integral of 
ipy could be evaluated when calculating the Franck Condon factor* 

Here n is the quantum number of Che level in the upper electronic state* 

Results for I 2 

The curves in figure 2 are for with v ■ 0 and n ■ 45* The wave 
functions can be checked by counting the number of crossings of the axis 
and 45)* Timnellng at the limits is apparent, and plots of showed 

that Che probability of finding the atom at the turn around points li the 
upper level is about four times greater than at the center* 

Although quantum numbers as high as 45 could be used for the upper 
level, for I 2 only v ^ 10 could be used for the lower level, otherwise 
too many crossing occurred* For v * 20 and 50, there were 21 and 57 
crossings, respectively* 

A second source of error occured for r < R. where the computer 
clamped the to zero if it diverged* In figure 2(b), the wave 
functions do not approach zero asymptotically as r ~ in this section. 
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N«w Afympcoclc Solution 


The failure of the wave functiona to approach zero asymptotically for r 
♦ - occurs in the region where the Franck Condon factors are calculated* 
This region is susceptible to the accumulated errors of the reiterations* 
Attempts to use fewer steps were unsuccessful* Accordingly, an asymptotic 
solution for i|>(r - “ ) was sought, with the view of carrying out a similar 

reiterative process with dr positive* 

For r « < r^ , the eiq^onentia' term in equation (3) dominates and 

the equation cakes the form 

y" - a^ exp (-2bx) y ■ 0 (13) 


where a,b are constants* By the transformation 


z 



(14) 


equation (13) becomes 



This is a form of Bessel's equation whose solution is 



(15) 


(lb) 
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wh«rt A And B AfA conttAntA dACtrainAd by t^A boundAry condlclonA, And 
^Of ^ *^A oodlfiAd BaaaaI functionA of chA fltAt and AACond kind* For 
largA z, 1q (A^/b) IncrAAAAA, and v OAAd y ^ 0; hancA, A ■ 0. Tha con- 
Atant B lA AvaluACAd hara by cha noraalizaclon procadura* For larga z, 
Kq haa an azymptoclc limic: 



(17) 


Tha anymptoclc aolutlon 1 a 



vhara A ■ S«9S1 x u Da and 

1 1 ^ j. 1 / \l 

^ ♦ ig. + y (ri - r)] 


(18) 


(19) 


Using aquACion (18), tha inlclAl value of r was astloatad whan ij* ■ 
lO’^o and tha raltaratlon procass carried out as before in 2,490 steps 
with dr positive* Tha plots of (fig* 3a) are for v ■ U, 10, and 45* 

A comparison of tha new method with tha previous shows tha two methods yield 
almost identical results batwaan tha classical turning points,' but tha now 
method gives values of ^ which approach zero asymptotically for r < B.* 

For r > K 4 ., thw new method gives diverging values of as tha program 
monitor was not applied hare to make » 0* With v - 10, tha error In 
near would only contribute to the normalization constant* However 
both methods gave a wave function which crossed the axis 50 times when v ■ 
45. 



Errors Using Hlghsr Quancuo labors 

Ic can bs stsn that for cha lowar lava I, both asynpcotlc solutions 
yiald functions which cross cha axis coo many cimas (fig* 3a) » buc cha func- 
cions bahava proparly for cha uppar laval (fig* 2a) up Co n ■ 43* In chasa 
plocs» valuaa of and w«y«' wara cakan co ba saro» as chay 

wara hoc Includad in rafaranca 2* In cha casa ot cha uppar laval, was 
calculatad from aquacion (S), which assumas (i»«y<^ * u* On calculacing 
from w^/4 b>cx« for Cha lowar lava! of l 2 , cha ascimatad valua 
was 2.32 aV, wharaas cha givan valua waa 1*54 aV. Using cha valua • 

2.32 aV for cha lowar laval, cha wava funcciop for v • 5u crossad cha axis 
50 cimas (fig* 4). Evidancly cha valua of w^y^ is important for high 
quantum numbars ond it cannot ba aasumad zaro* 

A valua of - d. 95 K iJ* cakan from rafaranca 1 was triad, 

buc again for v > 10, cha wava functions misbahaved* Attampcs wara mada co 
ascimaca o cha liaicing value of £(v) from aquacion (4) for 

larger v should ba aqual to If ■ 0, than cha limiting 

valua is D^' - w^/4 m«y«. Thus, fron aquacion (4), D« * I>a' “ 
m«y^ V|]^ whara vq is Cha quantum numbar Just whan dissociation 
occurs, which is givan (raf* 1) by 


v» ■ m 
D 2 tat X 

a a 


Uanca 


tat 


e^a 




( 20 ) 


( 21 ) 


Using the estimated valua tat^y^ ■ -0.0011b for v ■ 15, the wave function 
crossed the axis lb cimas, so the valua is not sufficiently accurate* 



Nert was estimated at the turning point for the level. For 

a given v, E(v) was calculated assuming (Dgye ■ 0, and the value of v 
at the turning point A (fig. 1) was calculated. The dotted curve is 
obtained assuming Dg* ■ '*‘e^e* length AB should be the con- 

tribution of the cubic term in equation (4) and, for v ■ 50, AB was 
approximately 2,000 cm“^; hence w^y^ ■ 0.01553. The value of E(v) 
was then recalculated, but the level was too low, and the wave function 
crossed the axis only 41 times. Hence this method also failed. Further 
study to estimate o^gFe ^ point on the Morse curves corresponding to 
the exact quantum level in question is underway. 

Orthogonality 

One check of the accuracy of the wave functions is to check for orthog- 
onality by evaluating the quantity t: 

e » ♦ dr 0 (22) 

N,v N,v' 

where v,v' are different vibrational quantun numbers for the sasie elec- 
tronic state. The wave functions are plotted for I 2 for 0 < v 30 in 
figure 5. The integration range was fixed from r • 3.318 to r ■ 2.300 A. 

The upper limit of r was determined when ^ was just greater than 10~ 
when v a 0. Even though, when v * 30, iteration started at a value of r 
where the asymptotic solution was doubtful, nevertheless the wave function 
appeared reasonable. The quantity e which 'Should be zero is shown in 
table 1. The values ranged from about 4 * 10*^ for the 0 to 1 case, but as 
v' became larger, e became larger and, for the v ■ 0 to 10 case, e was 
6 10 For adjacent levels, v to v ♦ 1, e < 10”^ up to v ■ 10. The 

low values of e imply that wave functions for 0 < v < 10 should be 
trustworthy in the lower level of I 2 . 


A 

■ 


3 


f 


>1 


10 




Franck Condon Factors 


Once the wave functions are known, the Franck Condon Factors, F, were 
obtained from 



(23) 


with the limits of Integration A,B and A,C chosen where the wave func'* 
tlons were nellglbly small* Ihe integration was carried out over 249 points 
using Slmspon's rule* Values for the transitions from v ■ 0, 1, *** 10 to n 

■ 0 to 25 were obtained for Nbe* The values for the v ■ 0 level are given 
In table 2* The energy difference between the lower and upper levels was 
expressed In terms of the wavelength of the absorbed photon, and hence plots 
of F vs* A were obtained* 

The vibrational levels in the lower electronic level are fairly widely 
spaced, and the relative populations obey Boltzmann Statistics Ny/No 

■ exp [-(Ey - Eo)/kT. The fractions Nv/Nq denoted by Qj for T ■ 300 

K and 0; for T * 1,000 K are given In table 3* Multiplying by Qi and Oj is 
equivalent to giving a statistical weight to the contribution of the various 
V levels in any absorption process* The Franck Condon factors luultiplied 
by 0^ and for NSe are plotted as if they were a continuous function of 
wavelength A in figures 6(a) and (b)* The peak values of F are about 
the same for v * 0 to 10, but on multiplying by ,Q 2 it can be seen that 
the contributions from the lower vibrational levels are the most 
important* 

In the case of I 2 , It has been shown that using the calculated rather 
than the measured value of gives different values to E and F 

(tables 1 and 4). The values of F for the two cases are compared for v 

■ 0 in figure 7 and the differences are small* The differences would be 
greater for higher v* In figures 8 and 9, plots are shown for I 2 for v ■ 0 
to 5 to the n ■ 0 to 45 levels assuming Dg ■ 2* 322b eV* The lower vibra- 
tional levels are much more closely spaced and the broadening at 1,UU0 K Is 
more evident* 


1 1 


Plots proportional to the absorption cross section per unit wavelength 
a a can be obtained from figures 6, 8, and 9 by oultlplylng F by the 
f requency v: 


o « vF 
a 


(24) 


where v Is the frequency of the absorbed photon and figures b» 8, and 9 
only roughly represent the absorption cross sectlons 4 

It should be emphasised that the selection rules have not been consld* 
ered here. However, It can be seen that at 1,000 K cUe contributions trout 
Che higher v levels become more Important; hence, the overall cross section, 
obtained by summing the Individual contributions, shows a broadening, as 
observed experimentally for diatomic molecules. A more accurate calculation 
would result If F were multiplied Instead by 

expj- (E^ - E^)/kx] 

J [exp(-E^ - E^)/kT] 
v-0 

which is Che true fraction out of Che total population. In such a case Che 
curves should not only broaden, but the peak value should also fall. 

Vibrational Sum Rule 

If the vibrational eigenfunctions are properly normalized, it can be 
shown from elementary properties of orthogonal functions that the sum of the 
squares of the overlap integrals (Franck Condon factors) summed over all 
values of the vibrational quantum numbers of the upper or of the lower 
states is equal to one (ref. 1): 


n/t! K ’"v '•'n **'1^ “ ^ 

V n 

The values of the Franck Condon factors for I 2 are still appreciable at 
the llnilclng value of n ■ 45 (table 4). In the case of NSe, with v ■ 0 the 
range of n froa 0 to 25 Is sufficient for a summation, as contributions from 
higher n values would be small* The sun for NSe equals 0*858~~sllghtly 
less than one* 


CONCLUSIONS 

Molecular wave functions have been calculated for several oiaterlals, 
and the methods can be applied In general* A new numeric asymptotic method 
of calculating wave functions Is presented which gives greater accuracy in 
the region of interest for calculating Frank Condon factors; otherwise, the 
two methods agree almost exactly both In amplitude and wavelength. The 
wave functions for the vibrational levels, v,n (in the lower/upper elec- 
tronic states, respectively) can be checked as they should cross the axis v 
and n times, respectively. For both methods of numerical calculation the 
results are accurate for V < 10 and work for n of at least 50 for I 2 * The 
use of the orthogonality theorem gave a method of checking the accuracy of 
the wave functions* The reason for the failure in the lower electronic 
level Is lack of information of the anharmonlc coefficient w^yg, which 
was entered as zero* In the upper level the “dissociation” energj’ D 2 was 
calculated on the assumption “ 0, and the vibrational levels Ejj 

were exact* In contrast. In the lower electronic level, entering m^yg 
as zero causes the Ugye(v + 0.5)^ term to be neglected; and this becomes 
Important at higher values of v* 

In practice, estimates of the wave functions are adequate for the 
lower V levels because, at normal temperatures, high v levels are hardly 
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populated* Attempts to estimate better values of ueye were unsuccess* 
ful* 

Plots were made of the Franck Condon factors plotted versus the wave* 
length of the absorbed photon, and the relative contribution from each of 
the lower v levels weighed with the population density* The resulting 
plots give an approximate representation of the absorption cross section per 
unit wavelength versus wavelength* Higher temperatures fill the upper v 
levels, resulting in broadening the absorption cross section, as observed 
experimentally* A check on the accuracy of the Franck Condon factors using 
the vibrational sum rule implied that the results were probably correct to 
within 20 percent for NSe. 


FUTURE WORK 


The present program, which uses an asymptotic solution at r > 
requires further work* Attempts to improve the accuracy of the wave func* 
tions at higher values of v should continue, by obtaining better estimates 
ui^ye and Cy* 

The method of solution, using the asymptotic value of at r < R. 

(see "New Asymptotic Solution”) should be further investigated* A minor 
modification is to monitor the values of <1* at r > R^., and if divergence 
occurs as r increases, then should be made either to decay exponen- 
tially (as it should) or else be clamped to zero* 

The new mthod has a number of advantages* It should give more accurate 
values of i|» at r < R., the region of interest in the upper level for 
calculating Franck Condon factors* In addition, the wave functions at 
values of Ejj < T« -t* (fig* 1) should be investigated* These colli- 
sions would result in dissociation with one of the atoms in an excited 
state, the events of greatest interest in the solar laser program* The 
problem of normalization is now different, as the wave function is a 


continuous wave as r increases to infinity* The possibility of “box” 
normalization arises (ref* 7)* 

In addition to dissociation probabilities for Horse curves* such as 
figure 1* where the upper level has a well* the possibility of calculating 
the transition to an upper level which is repulsive needs study* Repulsive 
curves can be represented by an exponential function asymptotically ap- 
proaching Dq. Hence* the method described under "New Asymptotic Solu- 
tion" would seem applicable* 

The Franck Condon factors should be applied to estimate absorption and 
emission coefficients* The expressions for the coefficients can be consid- 
ered to be in two parts* only one of which is given by the Franck Condon 
factor; the other part requires a knowledge of the dipole moments* Never- 
theless* the Franck Condon factors give an indication of the peak absorption 
frequency and bandwidths* 

The materials I 2 and NSe considered here were taken as examples of 
materials with closely and widely spaced vibrational levels* respectively* 
The programs should be applied to biatomic molecules which show promise for 
solar-pumped lasers* A list of possible candidates has been made alrea(^ 
based on estimates of solar* kinetic* and quantum efficiencies* 
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Table 1. Check of orthogonality for the lower electronic state of I 2 . (G 
should be zero.) 




D - 1.54238 ev 

p 

D ' - 2.3226 eV 

V 

v' 


6 

0 

1 

5. 31 3o3E-04 

4.52542E-04 

0 

2 

-8.59353E-04 

-1 .20373E-04 

0 

3 

-5.Sl3<»2E-04 

-1 . l23liE-03 

d 

4 

2.35<S7OE-03 

1 .3-9743E-04 

0 

? 

-2.32164E-03 

1 . 52555E-03 

0 

(» 

2. 31E47E-04 

-1 . 333O0E-03 

d 

? 

2.7SS70E-03 

-5. 13621E-04 

0 

3 

-5.09534E-03 

2. 30412E-0: 

0 


5. 55751E-03 

-2. 43934E-03 

0 

10 

-4. 3323oE-03 

1 . 01213E-O3 

e 

1 1 

1. 33l 3'?E-0:3 

1 . 3O457E-03 

0 

12 

3.30220E-03 

-3. 39540E-03 

0 

1 3 

-1.02070E-02 

4. 30O73E-03 

13 

14 

1.9O059E-02 

-3. 3r®>5E-03 

0 

15 

-2. 03337E-02 

2. 3O24.5E-03 

0 

16 

l.^*1093E-02 

2.111 15E-04 

0 

17 

-3. 7333OE-03 

-2. 93475E-03 

0 

13 

-4. 13057E-03 

5. 73323E-03 

0 

19 

2. 131 10E-02 

-3. 1 121 lE-03 

0 

20 

-4. 14329E-02 

1 . O2149E-02 

0 

21 

4. 29321 E-02 

-1 . 12771E-02 

0 


-4. 20040E-02 

1.2lolcE'’02 

0 

23 

1. 23o54E-02 

-1 . 23444E-02 


24 

-1 . 15173E-02 

1 . 2397-5E-02 

V 

0 

2 5 

3. 14-523E-02 

-1 . 2029 3E -02 


2-5 

-3. 9-5232E-02 

1 . 13303E-02 

0 

^ 1 

5. 15905E-02 

- 1 . IO72OE-02 

0 

23 

-7. 3o934E-02 

1 . 05751E-02 


29 

3 . 40150E-02 

-1.0O354E“Oi 

0 

30 

- 1 . 12730E-J31 

?. 379 74E-03 


10 

1 1 

1.331 39E-03 

1 . 3'94^7E-03 

10 

12 

3. 3022OE-03 

-3. 3354OE-03 

10 

13 

-1 . 0207OE-O2 

4. 30O7SE-03 

10 

14 

l.?^059E-02 

-3. 37695E-03 

10 

15 

-2.O3337E-02 

2. 3O24SE-03 
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Table 2. Values of Franck Condon factors assuming transitions from v - 0 
to the first 26 upper levels for NSe. 


LOWEP LEVEL 0 MO V= 0 


UFFEP LEVEL M 
0 
1 


•4 

e 

6 


10 
1 1 
12 

13 

14 

15 
1 E 

17 

18 
r? 
20 
21 
A. ^ 

23 

24 

25 


URVELENGTH-'R; ' 
4. 1 3205E^03 
4. O2436E+03 
3. 92388E>03 
3. 32994E-t-03 
3. 74198E+03 
3. 65947E+03 
3. 58197E>03 
3. 50908E 03 
3. 44043E+03 
3. 37569E+03 
3. 31459E+-03 
3. 25635E>03 
3.20225E+O3 
3. 15057E+03 
3. 10lo2E-*-03 
3.05521E>03 
3. 01 120E+03 
2. 96943E+03 
2. 92977E+03 
2. 89210E>O3 
2. 85631E+03 
2. 82229E-*-03 
2 . 7 w 3 9 4 E 0 3 
2. 75919E>03 
2. 72995E+03 
2. 702l5E-^03 


FRHHCL-CuhDON FRCTQP 
2. 43943E-03 
1 . 10533E-02 
2. 7681 4E-02 
4. 92627E-02 
7. 03649E-02 
3. 59509E-02 
9. 34244E-02 
■5. 2S4 40E-02 
8 . 59406E-02 
7. 52319E-02 
6. 29275E-02 
5.069 36E-02 
3. 96352E-02 
3. 02724E-O2 
2. 26549E-02 
1 . 66653E-02 
1 . 21274E-02 
8. 75237E-03 
6. 34492E-03 
4. 60472E-O3 
3. 24604E-O3 
2.21 126E-03 
1 . 53926E-0 3 
1 . 10463E-03 
7. 55805E-O4 
5. 33192E-04 


18 


Table 3 


Values of Ny/No for the lower electronic level of NSe at temper- 
atures of 300 K — Qi and 1,000 K — Q 2 . 


V 

0 

1 


4 


C' 

I 

S 

•4 

10 


Ql 

1 . ooooeEfOO 
1 . 07S62E-02 
I . 16341E-04 
1 . 25438E-tiE 
I . 35353E-03 
1 . 45994E-10 
1 . 57471E- 12 
1 . 8985 IE- 1 4 
1 . 33204E-18 
1 . 9 7 8 0 7 E - 1 8 
2. 1 3142E-20 


02 

1 . OOO00E+OO 
2. 570O8E-OI 
8. 80521E-O2 
I . 89758E-02 
4. 382S8E-03 
1 . 12123E-03 
2. 88177E-04 
7. 4O832E-05 
1 . 90347E-05 
4 . 8920 3E -08 
I . 25728E-08 
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Table 4. Franck Condon factors and wavelengths of absorbed photon for I 2 

(b) D - 1.54238 eV. 
e 


assuming (a) ■ 2.3226 eV and 


Upper 
Level N 

e 

1 

2 

3 

4 

5 

6 


24 

27 

26 

27 

23 

23 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 


Wavelength (a) 

6.35907E^03 
6.30925E>03 
6.260dlE^03 
6.2l369E<»-03 
6 . 1 6736E^03 
6. 12327E^03 
6.O7989E^03 
6. O3767E*03 


(a) 

2.3226 eV 


e 


(b) 

- 1.54238 eV 


5. 46927E^03 
5. 44322E>03 
5. 41787E>03 
5. 39320E^03 
5. 36919E^03 
5. 345S3E»03 
5. 323l0E>*>03 
5. 3O100E>O3 
5.27950£^O3 
5. 2536OE>03 
5. 23829E^03 
5.21355£^03 
5. 19937E»03 
5. 18074E4.O3 
5. 16265E>03 
5. 145O9E^03 
5. 128O5E<^03 
5. 1 1 152E>03 
5. O955OE>03 
5.07997E^03 
5.O6492E>03 
5. 05O36E^03 


4.01147E-11 
3.33556E-10 
1 . 04896E '03 
7. 63623b-08 
3.3S313E-07 
1 . 54644E>06 
5. 00323E-O6 
1 . 3513 3E-05 


I . 05341E-02 
1 . 30639E-02 
1.264O5E-02 
1.52894E-02 
1.44727E-02 
I . 51307E-02 
1. 74?3iE-w2 
1.59885E-02 
I.S0716E-02 
1 . 31736E-02 
I . 81219E-02 
I . 79234E-02 
I . 76334E-02 
1 . 54734E-02 
1 . 49269E-02 
1 . 45221E-02 
1 . 39S50E-02 
1 . 3 34O7E-02 
1 . 40 92OE-02 
1 . 5-‘-O?4E-02 
1 . l4i54E-02 
1 . 07S47E-02 


6.613I1E-11 
1.37 305E -O'* 
1.5044.jE-0i 
1 . 04427E-07 
5.08900E-07 
1 . 9807 6E-06 
6. 31746E” 06 
1 . 7 3429E-05 


s 

5.99659E-*-03 

3. 5O995E-05 

3.68480E-05 

9 

5.9566OE-*-03 

6. 33936E-05 

8.55731E-05 

10 

5.91763E^03 

1.23567E-04 

1 . 66596E-04 

1 1 

5.37979E^03 

2.42982E-04 

2. 72807E-04 

12 

5.34290E^03 

4. 05475E-04 

T. 1 3571E-04 

13 

5.30699E^03 

6.63923E-04 

27440E-04 

14 

5. 77203F->03 

1.01771E-03 

1 . 2 1 S ;■ *s E - 0 i 

15 

5. 73799E»03 

1 . 677O0E-03 

1 . 7 4 •; E - 0 i 

1 6 

5. 7O485E+03 

2. 40271E-O3 

2. 56 lE-0 •: 

17 

5.67257E>03 

2. 79347E-03 

. i 6 5 E - 0 

13 

5.64115E*03 

4. 23408E-03 

2 1 4 c £ - y 8 

19 

5.61 056E'^O3 

5. 261O0E-03 

Z . 504 30E-0 3 

20 

5. 53077E*03 

6. 52957E-03 

. y > 5 2 E - 0 i 

21 

5.55176E+03 

6 . SS834E-03 

. 0 1 1 >:■ 5 E - 0 i 

22 

5. 52353E^03 

8. 17924E-03 

9 . 4 4 ;j 0 2 E — 0 i 

23 

5. 49eO4E^03 

9. 35455E-03 

1 . 0 1 5 E - 0 4i 

i - ~ r - ^ 


1 . 20 3r?E-‘'-. 
1 . 2 ■? 1 8 7 E - 0 . 

1 . 6 3759E-02 
I . 706 1 IE -02 
1 . 7654IE-02 
1 . 7 9 7 6 7 E - 0 2 
1 .8227IE-02 
1 . 8274 iE-02 

1 . 6 2 ■? 7 6 E - 0 2 
1 .61512E-02 
1 . 5 ?•? 19E-02 
1 . 5'1 15E-02 
1 . 66757E-02 
: . 6 I 7 4 0 E - 0 2 
I . 546 I lE-02 

: . : 111 7e-o: 

1.2665IE-02 
1 . : 3326E-02 

• I » I 4 • * - /t * 

1 . 0 7 4 I 6 £ - 0 2 


0.376354 


0.3875S« 
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m 





r(A) 

(a) 



Figure 3. (a). Molecular wave functione for the lower electronic etate 

I 2 using a new asymptotic method; (a) v ■ 0, 10, 45. (b). Com- 
parison with previous asymptotic solution for v - 45 (continuous 
curve ■ new method: dotted curve ■ previous method). 




Figure 5 


8000 



R(flrJGST) 


Test of orthogonality for the first 30 vibrational levels in the 
lower electronic level of I 2 . D = lj^/4 uj x - 2 3226 eV 

q e g ee 
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50*Ql*F 


Figure 8 



Plot of 50 QF vs. wavelength for I 2 (v ■ 0 to 5; n - 0 to 45; lower 
level X* r*"; upper level H + u) at 300 K. 




Figure 9. Plot of 40 QF vs. wavelength for I 2 (v = 0 to 5 ; n * 0 to 45; lower 
level .X' r*” ; upper level n + u) at 1,000 K. 



